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Abstract 



Q^' A prescription for the short-distance expansion of Euchdean current cor- 

^ ■ relators based on a delocahzed modification of the multipole expansion of 

perturbative short-distance coefficient functions is proposed that appreciates 
the presence of nonlocal physics in the nonperturbative QCD vacuum. This 
^ ■ expansion converges better than the local Wilson OPE, which is recovered in 

the limit of infinite resolution. As a consequence, the usual local condensates 
in the Wilson OPE become condensates that depend on a resolution param- 
eter and that can be expressed as an infinite series of local condensates with 
increasing dimension. In a calculation of the nonperturbative correction to 
the ground state energy level of heavy quarkonia the improved convergence 
properties of the delocalized expansion are demonstrated. Phenomenological 
evidence is gathered that the gluon condensate, often being the leading non- 
perturbative parameter of the Wilson OPE, is indeed a function of resolution. 
The delocalized expansion is applied to derive a leading order scaling relation 
for /d//b in the heavy mass expansion. 



1 Introduction 



Establishing a connection between universal parameters of the nonperturbative QCD 
vacuum in Euclidean space-time on the one hand and an integral over weighted 
hadronic cross sections on the other hand has proven to be a very fruitful idea. [I] 
The three underlying assumptions of this approach are the analyticity of correla- 
tors of gauge invariant currents, the existence of a meaningful (asymptotic) large- 
momentum expansion of these correlators in the Euclidean region and the possibility 
of a dual description for hadron dynamics in terms of the dynamics of quarks and 
gluons. One may wonder whether these assumptions really are all independent. 
Given our poor analytic understanding of the nonperturbative quantum dynamics 
of quarks and gluons at present, we are not in a position to prove the analyticity 
assumption in QCD from first principles. On the other hand, one may, on general 
grounds and with the help of results from lattice gauge theory, try to connect the 
issue of a meaningful large-momentum expansion with that of quark-hadron duality. 

The conventional approach to the large-momentum expansion of the correlator of 
a gauge invariant QCD current j{x) is Wilson's operator product expansion (OPE)^, 

I / (i^xe^^- (Tj(x)jt(O) ) = E E ^Nim {O^m , (Q2 = _g2 > 0) . (1) 

J N=0 1=1 

In Eq. the contribution with = is solely due to perturbation theory, whereas 
> runs over the dimension of the local and gauge invariant composite opera- 
tors Om{0), and / labels operators of the same dimension. Operators ONih{0) and 
OATjiglO) with Ni < N2 need not be unrelated. Rather, ONihi^) be the result of 
gauge-covariant differentiations of the gauge invariant fc-point function correspond- 
ing to the content of fundamental fields in the local operator ONih{0). jT] In this 
sense On212{^) is a reducible operator. It is the purpose of this paper to exploit the 
relation of a chain of reducible local operators with nonlocal vacuum averages for 
constructing a modified version of the OPE. 

According to the standard interpretation jj, the physics of low-momentum, non- 
perturbative vacuum fluctuations resides in the vacuum averages (0^(0)), whereas 
high-momentum, perturbative propagation is contained in the Wilson coefficients 
CNiiQ'^)- If {Oni{0)) ~ then, apart from logarithmic factors arising from anoma- 
lous, perturbative operator scaling, we have cni{Q'^) ~ Q~'^- Here it is assumed 
that Q is the only relevant short- distance scale. Nonperturbative corrections to the 
= contribution are counted in powers of A/Q and are small if A <^ Q, and 
if In is sufficiently small. It is expected that the OPE is asymptotic at best, i.e. 
convergence is only apparent up to some critical value of N. 

^For the purpose of this introduction it is not necessary to consider the Lorentz and flavor 
structure of the current. Also, we assume that j{x) has no anomalous dimension. In general, the 
separation between high and low momenta in the OPE is ambiguous, and one has to introduce 
a factorization scale fi which cjv; and Oni depend upon. This dependence cancels in the whole 
series. 
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Figure 1: Schematic drawing of the short-distance function f{x) and the 2-point- 
correlator g{x) illustrating the scale hierarchy Af/Ag ~ A/Q ^ 1. 



The expansion into vacuum averages of local operators assumes correlations in 
the vacuum to be rigid, or in other words, it is assumed that the lower bound on the 
correlation lengths of the corresponding nonlocal gauge invariant fc-point functions 
is considerably larger than the short-distance scale Q^^. More specifically, one likes 
to think of correlation lengths being of the order ~ Aqqq, where Aqcd is the 
typical hadronization scale, i.e. it is commonly assumed that there is no hierarchy 
in nonperturbative scales. 

Given the exact knowledge of all gauge invariant QCD vacuum correlators and 
maintaining the possibility to factorize perturbative propagation from nonperturba- 
tive background fluctuations for each operator structure, one may view the OPE ^ 
as a degenerate multipole expansion of the perturbative part. To see this, we consider 
the simplified case of a nonperturbative and nonlocal dimension-4 structure corre- 
sponding to a slowly varying 2-point correlator g{x) falling off at distance ~ A~^ 
in its Euclidean space-time argument x. It is straightforward to generalize the sub- 
sequent considerations to arbitrary n-point correlators. In addition, we consider 
a short-distance function f{x), calculable in perturbation theory, that probes the 
vacuum at distance x ~ Aj ~ Q~^. Since the function f{x) describes short-distance 
physics at the scale Q ^ A it is strongly peaked compared to g{x) (see Fig.[T)). For 
demonstration purposes it suffices to consider an effectively one-dimensional prob- 
lem in which the chain of local power corrections to the purely perturbative result 
reads 



dxf{x)g{x) = dxdy f{y)5{x-y)g{x) 



dxdyfiy) 

L n=0 



2 



E 

n=0 



dyfiy) f 



f— 1 Y 
dx ^-^5^''\x) g{x) 
n\ 



f{k) 



E 

n=0 



dk 



k=0 



dk (ikr 



(2) 



The functions / and g are the Fourier transforms of / and respectively. From 
Eq. (j2D it follows that 



E 

n=0 



'-D 



X] 



(3) 



which is an expansion into degenerated multipoles. In electrodynamics the 3- 
dimensional generalization of Eq. (jS)) is used to determine the potential of a localized 
and static charge distribution p{x) at distances much larger than its spatial exten- 
sion (see e.g. Ref. pj). The first term on the RHS of Eq. (jH} corresponds to the total 
charge times the delta function located at the center of the charge distribution, the 
second term corresponds to the dipole moment times the derivative of the delta 
function, etc.. In the chain of local power corrections the derivatives of 5-functions 
generate higher- dimensional operators that are obtained by covariant differentia- 
tions of the correlator g{x) at the origin x = 0. We note that without specifying a 
regularization scheme the terms in the series of Eq. (j2I) do only exist, if the functions 
/ and g fall off sufficiently fast in the infrared and the ultraviolet, respectively. For 
g this statement translates to the property that the local operators, which are con- 
tained in the corresponding chain, do not have anomalous dimensions. The presence 
of anomalous dimensions implies a regularization scheme, which we will ignore in 
this introduction and in Sees. I2HS1 

Briefly switching to 4 dimensions, we now give an example for gauge invariant 
two-point correlations leading to the function g{x). The field strength correlator [S] 



where 



g^lVK\{x) = Tt i^g Gi,^{x) S{x,0)G^x{0) S{0,x) 
S{x,y) = pe^aJo'i-^M-) 



(4) 



(5) 



generates the following chain of vacuum expectation values (VEV's) involving local 
operators of increasing dimension^ 



Tr 6-^.(0) G'.a(O) 



Tr (g'G^MDpG^xm, Tr (g' G^,{0) DrDpG^x{0) 



(6) 
(7) 



This can most easily be seen in Schwinger's fixed point gauge x^y4^(x) = jj] if 
one assumes that the integration path in Eq. ((S)) is a straight line. In the limit, 

^The local operators with covariant derivatives are often expressed in terms of undifferentiated 
operators by appealing to the equations of motion and Biancchi identities. 
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where R = Ag/Af — >■ oo, the lowest power correction dominates, and operators 
such as in Eq. ((7)) are irrelevant in the OPE. The correlator of Eq. (0} depends on 
the closed path containing the points and x. For the purpose of lattice evaluations 
a straight line is always assumed. This may be adequate for the extraction of 
the highest mass nonperturbative scale being present in the corresponding gauge 
invariant correlations. 

The first VEV in Eq. is called the gluon condensate and VEV's with odd 
numbers of covariant derivatives do not contribute due to the experimentally well 
justified parity and time-reversal invariance of the QCD Lagrangian. Perturbatively, 
one can define the gluon condensate in such a way that its anomalous dimension 
vanishes to all orders in a^. [S| For the channels, which we will consider later, the 
anomalous dimension of the dimension six condensates are small. PQ We will there- 
fore ignore the problem of anomalous dimensions in these applications. We will 
come back to this issue in Sec. IHl 

Let us carry on assuming a one-dimensional world. Retaining only the lowest 
power correction at i? < oo corresponds to a saddle point approximation of the 
x-integral in Q. Expressed somewhat loosely in mathematical terms, the strongly 
peaked function / and the slowly varying function g lie in a dual space defined with 
respect to the bilinear form 

{f,g) ^ I dxf{x)g{x). (8) 

Formally, the set of functions 

e„(x) = ^^^^''\^) and e„(x) = , (9) 

which are used to construct the multipole expansion in Eq. (jH)), span a basis of the 
dual space, with the orthonormality relation 

The series of corrections to the saddle-point approximation, which involve the co- 
variant derivatives oi g{x) at a; = 0, is in general only asymptotic, and it is intuitively 
clear that the convergence properties worsen for R approaching unity. Much better 
convergence properties may be achieved by utilizing a multipole expansion of f{x) 
which starts out with a saddle point approximation based on an function of width 
Af instead of the infinitely narrow (5-function. In this paper we present an 

explicit construction of such an expansion and apply it to determine a modified ver- 
sion of the local OPE in Eq. (^. We refer to this expansion as "delocalized operator 
expansion" (DOE). 

The observation that the local expansion from Eq. breaks down for R ap- 
proaching unity and that the corrections to the saddle-point approximation can be 
substantial is well known and has led in the past to a number of phenomenological 
studies considering directly the non-local form in the first line of Eq. Q without any 
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expansion. 13 El E] In these studies the importance of including non-local effects 
has been demonstrated in a number of cases. In this approach the short- distance 
function f{x) has to be computed for arbitrary values of x and a specific ansatz for 
the form of the non-perturbative correlator g{x) has to be adopted. This leads to 
the unpleasant features that the determination of the perturbative short- distance 
function is considerably more complicated than for the Wilson coefficients of the 
local expansion, particularly a higher loop level, and that it necessarily involves 
some model-dependence which makes the estimates of theoretical errors difficult. 
The DOE has been constructed with the aim to provide an alternative formalism 
for non-local effects. The DOE can be applied with a definite power counting like 
the local OPE and keeps the computational advantages of the local expansion in the 
determination of the short- distance coefficients. The new feature of the DOE is that 
it relies on a generalization of the multipole expansion of the short- distance function 
f{x) based on functions having a width that can be adjusted to the short- distance 
scale Q. This leads to an improved convergence of the expansion for the price of 
having to introduce an additional parameter Q which we call "resolution scale" . The 
condensates in the DOE depend on Q and this dependence accounts for non-local 
effects. The DOE can be used to extract model-independent information on the 
non-local structure of the non-perturbative QCD vacuum with a well defined power 
counting (SecEJ, but also to simplify the determination of numerical predictions in 
a given model for the form of non-perturbative correlation functions g{x) (Secs.|3] 
and EI). 

The outline of the paper is as follows. In Sec. El we introduce the framework of 
the DOE assuming that the factorization of long- and short-distance fluctuations 
can be carried out without ambiguities. The delocalized expansion is associated 
with a resolution parameter Q and, for ~ Q, can lead to a series that has better 
convergence properties than the local OPE. Matrix elements in the DOE system- 
atically sum an inflnite series of local operators at each order in A/Q. At lowest 
order in A/Q the DOE implies an ^-dependent running gluon condensate. In Sec. El 
we review what is known about the gauge invariant gluonic fleld strength correlator 
from lattice simulations. In Sec.|3] the nonperturbative corrections to the ground 
state energy of a heavy quarkonium system are calculated in the local and the de- 
localized expansion in A/mv'^. Using a lattice-inspired model for the gluon fleld 
strength correlator, we show that the DOE has better asymptotic convergence than 
the local OPE. Sectional contains an extraction of the gluon condensate from ex- 
perimental data using charmonium sum rules and the Adler function in the V + A 
channel of light quark pair production. We flnd evidence that the gluon condensate 
is indeed running. We show that our results are consistent with the lattice mea- 
surements of the gluonic fleld strength correlator. In Sec. El we show how the DOE 
is constructed for the more general case that the factorization of long- and short- 
distance fluctuations is ambiguous and needs to be carried out in a regularization 
scheme. We then apply the DOE to derive a scaling relation for the ratio of meson 
decay constants /d//s at leading order in the heavy mass expansion. A summary 
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and an outlook are given in Sec.[7| Finally, in App.^we show that the gluon field 
strength correlator has approximately the same tensor structure as the local gluon 
condensate, which simplifies computations in the investigations in the main body 
of the paper. Appendix El contains analytical formulae for the local short- distance 
coefficients needed for the examination of quarkonium energy levels in Sec.|3] 



2 Delocalized operator expansion: framework, 
evolutions equations and power counting 

To construct a modified version of the multipole expansion used in Eq. (jH)) we need 
to find a generalization of the basis functions in Eq. (jU)) satisfying an orthonormality 
relation in analogy to Eq. (jlOj) . In Ref. JU] a lowest-order analysis based on a 4- 
dimensional spherical well was performed. Here, we construct a basis of the dual 
space by appealing to the orthonormality properties of the Hermite polynomials 
and their duals. For a treatment in Cartesian coordinates the following basis is well 
suited: 

nn+l 

\ rKTV. 



For even values of n the first three Hermite polynomials read 

Hoix) = 1, H2ix) =4x^-2, Hi{x) = 16x^-48x^ + 12. (12) 

The width of the basis functions is of order and one recovers the basis of 
the local expansion in the limit — > cxd: 

e„ — > Cn and e„ — ^ e„ . (13) 

An expansion of f{x)g{x) using the resolution dependent basis {ej^,ej^} generates 
the following representation of the x-integral in Eq. Q 

/oo 
dxf{x)g{x) = J2 fni^) 9n{n) ^ (14) 
-oo „ n 



n=0 



where 



(2^)" 
d 



{2ny 

6 



i/.(n(.^))/» 



(15) 



fc=0 



/on+1 
dx^-H^{nx)e-'''-'g{x 



^^e^9{k). (16) 

We call Q the "resolution parameter" . It is intuitively clear that the series in Eq. ()14|1 
has better convergence properties, if fl is of order Q, because in this case the term 
(/, 60)60(0;) accounts much better for the form of f{x) than for fl = 00. The terms 
/„(n) are the resolution-dependent short-distance coefficients, and the 5'„(fi) are the 
resolution- dependent condensates. Note that, if /(x) is positive definite within its 
main support < Q~^, then /n(oo) > for n not too large. On the other hand, 
if g{x) is positive definite and monotonic decreasing for x > 0, then gn{^) < Qni^') 
for VL < n'. 

The relation between basis functions for different resolution parameters Q and 
Q' can be obtained from the properties of the Hermite polynomials and reads 



m=0 

00 

e^ix) = Trunin', n)eZ{x), (17) 



m=0 



,^L^^ ( . „ „ ) ^ , — m > and even . 



where 



, otherwise . 

Equation (fTHj) leads to the following relations between the short-distance coefficients 
/„ and the corresponding condensates (?„ for different resolutions fl and fl': 



/„(fi) = ^ a„„(fi,fi')/m(f^') (19) 

m=0 

00 

gnin) = J2('--n{n',n)gmin'). (20) 
m=0 

The coefficients a„m(^^,fi') satisfy the composition property 

00 

Y am{^, fi') (^^m{^', = ^nU^, n") , (21) 

i=0 

and for equal resolution we have 

(^nm(fl: fl) ^nm ■ (^^) 
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Obviously, the transformations a„m(^)^') form a group. It is instructive to specify 
the relations in Eqs. (fT^ and (pUj) for Q' = oo, 

where the explicit expressions for fn{oo) and 5'„(oo) can be read off Eq. We 
see that each short-distance coefficient fn{^) can be expressed in terms of a finite 
linear combination of the local Wilson coefficients fi{oo) for i < n. In particular, 
the short- distance coefficient fo{^) and the Wilson coefficient of the first term in the 
local expansion coincide since Ho{Qx) = 1, i.e. the short- distance coefficient in the 
saddle-point approximation is f2-independent. On the other hand, the resolution- 
dependent condensates gni^) are related to an infinite sum of local condensates with 
additional covariant derivatives multiplied by f2-dependent coefficients. From now 
on we adopt the language saying that resolution-dependent condensates of a given di- 
mension can be expressed as an infinite sum of local condensates of equal and higher 
dimension. And similarly we say that each resolution-dependent short- distance co- 
efficient associated with a condensate of a certain dimension can be expressed as a 
finite sum of local short- distance coefficients associated with condensates of equal 
or smaller dimensions. 

Note that in Eq. (fT^ the term /o(cxD)(72(oo)/(2f2^) contained in /o(fi)(7o(^) is 
cancelled by the term — /o(oo)(72(oo)/(2f2^) in f2{^)g2{^) which exemplarily shows 
the reshuffling of the OPE. As we have argued above, the Wilson coefficients /n(oo) 
have equal sign for n not too large and therefore the short- distance coefficient f2{^) 
can be much smaller than f2{oo) for f2 ~ Q. This feature and the fact that gn{^) < 
Qnioo) under the conditions mentioned above can, at least for small n, severely 
suppress higher order contributions. The DOE of the current correlator in Eq. (^Q) 
has the same parametric counting in powers of A/Q as the OPE, as long as Q is 
not chosen parametrically smaller than Q. Taking into account that (?„ ~ A" and 
fn ~ <5~" in the local expansion, we also have 

9n{n) ~ A"^: (^)' ~ A", Un) ~ Q-^Ei^y - (24) 



/o(oo) 
/2(oo) 



/o(oo) 



3 3 



[n/2] 

E 

00 

E 



:n-2i)\i\ V 4^2 



/„_2i(oo) 



n + 2iy. 



gn+2i{oc) 



(23) 
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in the DOE as long as f2 3> Q- At finite resolution, however, one obtains additional 
summations of powers oi Q/VL and A/f2 in the short-distance coefficients and the 
condensates, respectively. Adopting a lattice-inspired model, we will see in Sec. IH 
that for r2 ~ Q the expansion in powers of A/Q has an additional suppression by 
powers of a small number as compared to the expansion ioi VL = oo. 

It is straightforward to extend the previous results to an arbitrary number of 
dimensions by applying the resolution-dependent expansion independently to each 
coordinate. This also accounts for the treatment of arbitrary ra-point correlations 
in an arbitrary number of dimensions. For convenience we consider here the same 
resolution parameter VL for each coordinate. It is straightforward to derive the 
resolution-dependent expansion for a more general situation. 

If / = /(xi, . . . , Xd) is a strongly peaked short- distance function in d dimensional 
Euclidean space-time the analogous expansion to Eq. reads 



d'^xf{Xi,...,Xd)g{Xl,...,Xd) = f m,...,na{^) gnu..., n^i^) , (25) 

ni,...,na=0 



where 



/?ii,...,nd(^) 



d'^xfiXi, ...,Xd) 



r d 

n 



Hn^nxi] 



n 



\ (2 ny^ V V dk 



9n,,...,n,{^) = d'^X \{^^HnXVLX^) 



f{ki, . . . , kd) 



e g{xi, 



(26) 



ki,...,ka=0 



d'^k 
(27r)« 



n 

i=l 



[ikn,] 



giki,...,kd) 



(27) 



The relation between the short-distance coefficients and the condensates in Eqs. ()26|) 
and (j^7j) for different resolution parameters is derived from Eqs. ()17p for each index 
Hi in analogy to the expressions in Eqs. ()19j) and (j20|) . As in the one- dimensional case 
each short- distance coefficient fni....,na{^) can be expressed in terms of a finite linear 
combination of the local Wilson coefficients fi^...,i^{oo) for ik < rik, and the short- 
distance coefficient /o,...,o(^) is resolution-independent. The resolution-dependent 
condensates (7„^^...^„^(f2) are related to an infinite sum of local condensates (/^^^...^^^(oo) 



for ik 

{xl + . 



> Uk- For example, assuming that (7 is a function of the distance 
. . + a;^)"*^/^ only, one finds 



\x\ 



5'o,...,o(^) 



2Q2 



5'2,...,o(ooj + • • • +5'0,...,2(OO) + 



-,dpd,g{Q) 



{2i 
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It is an easy exercise to derive evolution equations for the resolution- dependent 
short- distance coefficients and the condensates. In the one-dimensional case the 
evolution equations are obtained directly from Eqs. (fT^ and (^0)) and read 

^/"(") - ^ 

d {n + l){n + 2) 
j^9n{n) = g„M^). (29) 

Equations ()29p imply that there is only mixing between terms that differ by two 
mass dimensions. This is a consequence of the specific choice of the basis functions 
in Eqs. (jllj) . which are constructed from the Hermite polynomials and the Gaussian 
function. For the short- distance coefficients lower dimensional terms always mix into 
higher dimensional ones. Prescribing short-distance coefficients at i7 = oo, the solu- 
tions are given by Eqs. ()23|) . For the condensates higher dimensional terms always 
mix into lower dimensional ones. Their evolution equations are much more interest- 
ing because they show that knowing the leading condensate go{^) implies that all 
subleading condensates of higher dimensions are obtained by differentiations with 
respect to Q. The analogous features also exist for the d-dimensional generalization 
of the evolution equations, which have the form 

d _ {m + i){m + 2) 

K + l)K + 2) 

• • • fl'ni,...,nd+2l"j • [OO) 



3 The gluonic field strength correlator 

Since it is the most important 2-point function, which hence we will heavily draw 
upon in subsequent sections, we review here what is known from lattice calculations 
about the gauge invariant, gluonic field strength correlator (jl)). Without constrain- 
ing generality the following parametrization of this function in Euclidean space-time 
was introduced in Ref. 

To separate perturbative from nonperturbative contributions the scalar functions D 
and Di are usually ffited as 
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Z,.(,=) = ^.exp(-M) + ^e.p(-M,. (32) 



The power-like behavior in Eqs. ()32j) at small |x| is believed to catch most of the 
perturbative physics, although it is known that partially summed perturbation the- 
ory may generate ultraviolet-finite contributions for = as well. However, 
magnitude estimates are difficult because the integration over the corresponding 
poles in the Borel plane is known to be ambiguous. [TTl We ignore this subtlety 
and assume that the purely exponential terms in Eqs. ()32j) exclusively carry non- 
perturbative information. In a recent unquenched lattice simulation ^2] (see also 
Ref. ^3]) where the gluon field strength correlator was measured with a resolution 
of ^ (0.1 fm)~^ ^ 2 GeV between 3 and 8 lattice spacings it was found that 

^ ^ 9 and ^ 0.7 GeV. (33) 

It is conspicuous that the inverse correlation length is somewhat larger than the 
typical hadronization scale Aqcd- Whereas the actual size of Aq and Ai depend 
quite strongly on whether quenched or unquenched simulations are carried out and 
which values for the light quark masses were assumed, the ratio Aq/Ai and the 
correlation length Xa were found to be quite stable. [12] 

Because Aq ^ Ai we may neglect the Di contribution in Eq. ()31|) at the level of 
precision intended in our subsequent analyses and write 

V^a''"'*(^) = ^ ( '^M'^^^A - 5f,x5uK ) g{\x\) . (34) 

This simplifies the computations involved in the DOE since in this approximation 
the gluon field strength correlator has the same x-independent tensor structure as 
the local gluon condensate 

,2 



vr 



d^^"(o) = ( s,j., - s,,s.^ ) Y (0 ^ g;mg;m o) . (35) 



a. 



TT 



It is shown in App.^that the lattice-implied dominance of the tensor structure in 
Eq. ()34|) is consistent with a determination of 



TT{g'G,MDaDpG^x{0)) (36) 

using equations of motion, Biancchi identities, and phenomenologically obtained 
values for the VEV's (g'jX) and (gfabcGlAGj ^ 



We note that due to the cusp of the exponential ansatz in Eq. , higher deriva- 
tives at a; = have nonlogarithmic UV-singularities that can only be defined in a 
regularization scheme. So assuming an exponential behavior down to arbitrary small 
distances would cause inconsistencies in the Wilson OPE for operators containing 
covariant derivatives, which are known to have logarithmic anomalous dimensions. 
Since unquenched lattice simulations have not been carried out for distances below 
3 lattice spacings, the simple exponential ansatz is invalid for some region around 
X = with Ixl < Xa- 
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4 Nonperturbative corrections to the heavy 
quarkonium ground state level 

Among the early applications of the OPE in QCD was the analysis of nonperturba- 
tive effects in heavy quarkonium systems. |141 ITH] Heavy quarkonium systems are 
nonrelativistic quark-antiquark bound states for which there is the following hierar- 
chy of the relevant physical scales m (heavy quark mass), mv (relative momentum), 
mv'^ (kinetic energy) and Aqcd: 

m > mv > mv"^ > Aqcd • (37) 

Thus the spatial size of the quarkonium system ~ {mv)~^ is much smaller than the 
typical dynamical time scale ~ (mw^)"^. We note that in practice the last of the con- 
ditions in Eq. ()37|). which relates the vacuum correlation length with the quarkonium 
energy scale, is probably not satisfied for any known quarkonium state, not even for 
T mesons Only for top-antitop quark threshold production condition (|H7j) is 
believed to be a viable assumption. [T7j 

In this section we demonstrate the DOE for the nonperturbative corrections to 
the n'^'^'^^Lj = l^Si ground state for different values of m in the model defined 
in Eq. ()42|1 . We adopt the local version of the multipole expansion (OPE) for the 
expansion in the ratios of the scales m, mv and mv"^. The resolution dependent 
expansion (DOE) is applied with respect to the ratio of the scales mf^ and Aqcd- 
The former expansion amounts to the usual treatment of the dominant perturbative 
dynamics by means of a nonrelativistic two-body Schrodinger equation. The interac- 
tion with the nonperturbative vacuum is accounted for by two insertions of the local 
xE dipole operator, E being the chromoelectric field. J3] The chain of VEV's of the 
two gluon operator with increasing numbers of covariant derivatives times powers of 
quark-antiquark octet propagators i.e. the expansion in A/mv"^, is treated in 
the DOE. Our examination is not intended to represent a phenomenological study 
of non-perturbative effects in heavy quarkonium energy levels, but to demonstrate 
the result of the DOE for a natural choice of the resolution scale in comparison to 
the local expansion and the exact result in a given model. 

At leading order in the local multipole expansion with respect to the scales m, 
mv, and mv"^ the expression for the nonperturbative corrections to the ground state 
energy in a form analogous to Eqs. © and (fT^ reads 

/oo 
dtfit)g{t), (38) 
-co 

where 

fit) = ^ / ^ / d'^ I d'y m (xy) Go (x, y, - go) m , (39) 

with 

GoUy,--) = f:(2/ + l)p/^)G/x,y,--), 
V m/ ^ \xy/ \ mj 
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Figure 2: The perturbative short-distance function f{t) for m = 5 GeV and as = 
0.39 (sohd hne) and the leading term in the delocahzed multipole expansion of f{t) 
ioT Q = k"^ /m (dashed hne). 



mj 271 ' ' ^ f^^{s + l + l-^){s + 2l + l)\ 



^3/2 



s=o y-^ ^ " ^ ^ 12 k 



k = -ma^ (40) 

The term Go is the quark-antiquark octet Green- function ^^1; and (j) denotes the 
ground state wave function. The functions P„ and are the Legendre and Laguerre 
polynomials, respectively. Since we neglect the spatial extension of the quarkonium 
system with respect to the interaction with the nonperturbative vacuum, the inser- 
tions of the xE operator probe only the temporal correlations in the vacuum, which 
effectively renders the problem one-dimensional. We note that t is the Euclidean 
time. This is the origin of the term it in the exponent appearing in the definition 
of the function /. For / it is implied that possible odd components are removed 
through the operation f{t) — > (/(t) + /(— 1))/2 because g is an even function. In 
Fig. 121 the function f(t) (solid line) is displayed for m = 5 GeV and = 0.39. Since 
the spatial extension of the quarkonium system is neglected and the average time 
between interactions with the vacuum is of the order of the inverse kinetic energy, 
the characteristic width of / is of order (mf^)~^ ~ (ma^)~^. As a comparison we 
have also displayed the function [/ (it'/(t')eQ (t')]eQ (t) for Q = k'^/m (dashed line), 
which is the leading term in the delocalized multipole expansion of /. The values of 
the first few local multipole moments fn{oo), which correspond to the local Wilson 



13 



coefficients, read 



/o(oo) = 1.6518 



36 k* ' 



m 



/2(oo) = 1.3130 




36 fcS 



/4(oo) = 7.7570 




/6(oo) = 130.492 




36fci6 



(41) 



/8(oo) = 4474.1 




/lo(oo) = 262709.3 




The term fo{oo) agrees with Ref. [1311^1 and f2{oo) with Ref. [TH]. The results for 
/n>2(oo) are new. Analytical expressions for the numbers shown in Eqs. 1)411) are 
given in App.lBl 

Let us compare the local expansion of -E"^ with the resolution-dependent ex- 
pansion using the basis functions of Eq. (|TT)l . For the nonperturbative gluonic field 
strength correlator we use a lattice inspired model of the form 



This model has an exponential large-time behavior according to Eq. ()32p and a 
smooth behavior for small t. The local dimension 4 gluon condensate in this model 
is 



We note that for the purpose of this examination the exact form of the model for the 
nonperturbative gluonic field strength correlator is not an important issue as long 
as the derivatives of g{t) at t = are well defined, see e.g. Refs. [UllTn] for different 
choices of models. 

In Tab.^ we have displayed the exact result in the model ()42p and the first 
four terms of the resolution-dependent expansion of E"'^ for the quark masses m = 
5,25,45,90, 175 GeV and for Q = oo and Q = k'^/m. For each value of the quark 
mass the strong coupling has been fixed by the relation = as{k). We note that 
the series are all asymptotic, i.e. they are not convergent for any resolution. The 
local expansion {Q = oo) is quite badly behaved for smaller quark masses because 
for k'^/m < A^^ any local expansion is meaningless. In particular, for m = 5 GeV 
the subleading dimension 6 term is already larger than the parametrically leading 
dimension 4 term. This is consistent with the size of the dimension 6 term based 
on a phenomenological estimate of the local dimension 6 condensate in Eq. ()36p . see 
App.|Xl For quark masses, where fc^/m > A^^, the local expansion is reasonably 
good. However, for finite resolution Q = k'^/m, the size of higher order terms is 
considerably smaller than in the local expansion for all quark masses, and the series 
appears to be much better behaved. The size of the order n term is suppressed by 
approximately a factor 2^" as compared to the order n term in the local expansion. 
We find explicitly that terms in the series with larger n decrease more quickly 
for finite resolution scale as compared to the local expansion. One also observes 



g{t) = 12Aoexp[~^t^ + X\/XA + l), 
Ao = 0.04 GeV\ A;^^ = 0.7 GeV, 



(42) 




(43) 
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0. 1 


o.o 










6 


—00.4 


Add 
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90 


0.17 


1.156 


1.05 





1.15 


1.15 


1.07 


1.07 










2 


-0.17 


0.98 


-0.04 


1.02 










4 


0.18 


1.16 


0.04 


1.07 










6 


-0.65 


0.51 


-0.07 


1.00 


175 


0.15 


1.750 


0.245 





0.258 


0.258 


0.249 


0.249 










2 


-0.016 


0.242 


-0.005 


0.244 










4 


0.008 


0.250 


0.003 


0.247 










6 


-0.012 


0.237 


-0.003 


0.244 



Table 1: Nonperturbative corrections to the heavy quarkonium ground state level 
at leading order in the multipole expansion with respect to the scales m, mv and 
mv^ for various quark masses m based on the model in Eq. ()42j). Displayed are the 
exact result and the first few orders in the DOE for f2 = oo and = k^ jm. The 
numbers are rounded off to units of 0.1, 0.01 or 0.001 MeV. 



that even in the case k^ jm < A^^, where the leading term of the local expansion 
overestimates the exact result, the leading term in the delocalized expansion for 
Q = k'^/m agrees with the exact result within a few percent. It is intuitively clear 
that this feature is a general property of the delocalized expansion, and we believe 
that it should apply to any problem for which the local expansion in the ratio of 
two scales breaks down because the ratio is not sufficiently small. 

We would like to emphasize again that the previous analysis is not intended 
to provide a phenomenological determination of nonperturbative corrections to the 
heavy quarkonium ground state energy level, but rather to demonstrate the DOE 
within a specific model. For a realistic treatment of the nonperturbative contribu- 
tions in the heavy quarkonium spectrum a model-independent analysis should be 
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carried out. In addition, also higher orders in the local multipole expansion with 
respect to the ratios of scales m, mv and mv"^ should be taken into account, which 
have been neglected here. These corrections might be substantial, in particular for 
smaller quark masses. However, having in mind an application to the bottomonium 
spectrum, we believe that our results for the expansion in A/mv'^ indicate that going 
beyond the leading term in the OPE for the bottomonium ground state is proba- 
bly meaningless and that calculations based on the DOE with a suitable choice of 
resolution should be more reliable. 



5 Extraction of the running gluon condensate 

The resolution-dependent condensates are either determined phenomenologically 
from experimental data or from lattice measurements. It is the purpose of this 
section to extract the resolution-dependent dimension 4 condensate from an ana- 
lysis of moment ratios for the charmonium system and a sum rule for the Adler 
function using the V + A spectral function measured at LEP in hadronic tau de- 
cays. In contrast to the demonstration in the previous section, where within our 
approximation, the vacuum was probed only by the temporal dynamics of the heavy 
quarkonium system, the charmonium moments and the V + A sum rules probe the 
full four-dimensional space-time structure of the vacuum. We use the 4-dimensional 
version of the DOE with a common resolution scale for each coordinate as given in 
Eqs. fl25|l - (j^H|l . For convenience we define 



and in the following we will refer to it as the "running gluon condensate". For 
Q = oo the running gluon condensate coincides with the local gluon condensate 



The determination of the local gluon condensate G'^j{oo) from charmonium sum 
rules was pioneered in Refs. [20] • By now there is a vast literature on computations 
of Wilson coefficients to various loop-orders and various dimensions of the power 
corrections in the correlator of two heavy quark currents including updated analyses 
of the corresponding sum rules (see e.g. Refs. |2IIE3)- Because the short- distance 
coefficient of the running gluon condensate is fi- independent, we can use the same 
strategy as in previous analyses, where the local condensate was determined. The 
only difference is that we have to keep track of its dependence on the characteristic 
short- distance scale. 

The relevant correlator is 



iq,q. - q^9,u) n^(Q') = i j d'xe''^^ {Tf,{x)jm) , Q' = -q' , (45) 




(44) 



{iaj7r)GU0)G;M)- 



5.1 Charmonium sum rules 
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where = C7^c, and the n-th moment is defined as 

Assuming analyticity of H'^ in away from the negative, real axis and employing 
the optical theorem, the nth moment can be expressed as a dispersion integral over 
the charm pair cross section in e^e~ annihilation, 

where s is the square of the cm. energy and Qc = 2/3. We consider the ratio 

Mn-1 



(48) 



and extract the running gluon condensate as a function of n from the equality of 
the theoretical ratio using Eq. ()46|) and the ratio based on Eq. ()47p determined from 
experimental data. 

The analytic form of the ra-th moment in the local OPE up to terms with di- 
mension 6 reads [21] 



^(6) 



(n + 3)! 1 



+ 



(n-l)!(2n + 5) 9 (4m2)2 ' 

2 (n + 4)! (3n2 + 8n - 5) 1 
45 (n- l)!(2n + 5) (2n + 7) 9 (4m2)3 ' 



(6) _ _ 8 (n + 2)!(n + 4)(3n3 + 47n2 + 244n + 405) 1 

~ ^135 (n-l)!(2n + 5) (2n + 7) 9 (4m2)3 ' ^ ^ 

where = (^?V'''^G'^,Gt,Gi^) and ((^^j^) ^ (^^j-j-), being the light fiavor 

singlet current. 

In contrast to the calculation of the quarkonium energy level in Sec.|31 where we 
had a nonrelativistic powercounting argument allowing for consistently considering 
only two gluons interacting with the vacuum, the local OPE of the charmonium 
moments is inconsistent if only the interaction of two gluons with the vacuum is 
accounted for. This is because the Wilson coefficients for condensates of dimension 
6 and higher are only gauge-invariant if the interaction of any number of gluons that 
can contribute at a certain dimension is taken into account. El Therefore, a reliable 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


p = l 


0.13 


0.36 


0.56 


0.73 


0.87 


0.99 


1.08 


1.16 


p = 2 


0.09 


0.25 


0.42 


0.57 


0.69 


0.80 


0.89 


0.97 



Table 2: Ratio of the local dimension-6 contributions contained in the running gluon 
condensate and in the full OPE for r„ as a function of n. 



extraction of the running gluon condensate is only possible if the higher dimensional 
local terms summed in 5^^\g'^G'^){VL) represent a reasonably good approximation to 
the full set of terms with the corresponding dimensions in the local OPE. The most 
important subleading contributions in the local OPE are the dimension 6 terms 
shown in Eq. (j49j) . Using relation (|65|) . one can derive that 

{g^GD^G) ^ {g'G,MDpDpG,um = 2{g'f) -2{g'fG'). (50) 

Assuming further Eq. (|34j) . then the local dimension 6 contribution contained in 
reads 

^ {g^GD'G) (51) 

In Tab. 121 we have displayed, for Q = 2mc/n, the ratio of the sum of local dimension 
6 contributions in the running gluon condensate correction according to Eq. ()51|) 
and the sum of the local dimension 6 term in Eq. ()49j) for r„. We have used the 
standard values {g^fG"^) = 0.045 GeV^ obtained from the instanton gas approx- 
imation P and {g^f) = -pA/3{AiT)'^al{qq)'^ with = 0.7 {p ^ 0.7 GeV) and 
{qq) = —(0.24 GeV)'^. The choice p = 1 corresponds to exact vacuum saturation. 
As explained later we believe that fl = 2mc/n is an appropriate choice for the reso- 
lution scale. Indeed, one finds for the relevant ratios r„ (2 < n < 7) that the local 
dimension 6 contribution contained in the running gluon condensate has the same 
sign and roughly the same size as the dimension 6 power correction in the full OPE 
for larger n, where higher dimensional contributions are more important. 

For the experimental moments we use the compilation presented in Ref. j23j . 
where the spectral function is split into contributions from the charmonium reso- 
nances, the charm threshold region, and the continuum. For the latter the authors 
of Ref. j2Sl used perturbation theory since no experimental data are available for the 
continuum region. We have assigned a 10% error for the spectral function in the 
continuum region. We believe that this is sufficiently conservative since a large part 
of the error drops out in the ratio r„. For the purely perturbative contribution of the 
theoretical moments we used the compilation of analytic 0(af ) results from Ref. j2Sl 
and adopted the MS mass definition mc(mc) (for any renormalization scale /x). For 
the Wilson coefficient of the gluon condensate we have used the expression given in 
Eqs. fl49|) . We have checked that for n < 8 the perturbative Ola"^) corrections do 
not exceed 50% of the 0{as) corrections for /i between 1 and 4 GeV. 
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Figure 3: The running gluon condensate as a function of n when extracted from 
the ratio of moments r„ for rndrfic) = 1.23 (white triangles), 1.24 (black stars), 1.25 
(white squares) and 1.26 GeV (black triangles). The area between the upper and 
lower symbols represents the uncertainties. The thick solid line indicates the running 
gluon condensate as it is obtained from the lattice-inspired ansatz in Eqs. ()32j) 



Our result for the running gluon condensate as a function of n is shown in 
Fig. El for n < 8 and mcirfic) = 1.23 (white triangles), 1.24 (black stars), 1.25 
(white squares) and 1.26 GeV (black triangles). For each value of rnd^c) and n 
the area between the upper and lower symbols represents the uncertainty due to the 
experimental errors of the moments and variations of as{Mz) between 0.116 and 
0.120 and of /x between 1 and 3 GeV. We see that the running gluon condensate 
appears to be a decreasing function of n. Let us compare this behavior of the running 
gluon condensate with the one obtained from the lattice-inspired exponential ansatz 
for g{x) in Eq. (jH^ using the definition of Eq. (fT^ .^ 

Here, erf is the error function. Due to the exponential behavior of the lattice-inspired 
ansatz for g{x) the expression in Eq. ()52|) is a monotonic increasing function of Q. 
Since the exact width of the short-distance function / for the moments is unknown 
we estimate it using physical arguments. For large n, i.e. in the nonrelativistic 
regime, the width is of the order of the quark cm. kinetic energy mv'^, which scales 

■^As explained in Sec. El Eq. I|52(l corresponds to a computation of the running gluon condensate 
in a regularization scheme with a cut-ofF of order Q. 
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like m/n because the average quark velocity in the n-th moment scales like Xj ^fn. 
For small n, on the other hand, the relevant short-distance scale is just the quark 
mass. Therefore, we take f2 = 2mc/n as the most appropriate choice of the resolution 
scale. In Fig. El the running gluon condensate from Eq. ()52j) is displayed for A^^ = 
0.7 GeV and, exemplarily, Aq = 0.04 (i.e. {{as /7r)G'^)i^tioo) = 0.024 GeV^) and for 
fl = (2.5 GeV)/n as the thick solid line. We find that the ra-dependence of the lattice 
inspired running gluon condensate and our fit result from the charmonium sum rules 
are consistent. This appears to confirm the lattice result for the correlation length 
A^. However, the uncertainties of our extraction are still quite large, particularly 
for n = 2, where the sensitivity to the error in the continuum region is enhanced. 
We also note that the values for the running gluon condensate tend to decrease with 
increasing charm quark mass. For mcirric) = 1.27 GeV we find that ((as/7r)G^) 
ranges from —0.01 GeV^ to +0.01 GeV^ and shifts further towards negative values 
for mc(mc) > 1.28 GeV for all n > 2. Assuming the reliability of the DOE as 
well as the local OPE in describing the nonperturbative effects in the charmonium 
sum rules and that the Euclidean scalar functions D{x) and -Di(x) in Eqs. ()31|) ) are 
positive definite, our result disfavors mcirfic) > 1.28 GeV. 

5.2 V + A sum rule 

Recently, the spectral function for light quark production in the V + A channel has 
been remeasured by the Aleph [2^] and Opal collaborations from hadronic r 
decays for < m^. In the OPE the associated current correlator is dominated 
by the gluon condensate, and the dimension 6 power corrections that are not due 
to a double covariant derivative in the local gluon condensate are suppressed. The 
corresponding currents are J^^^ = ^t7/x(l ± 75)*^ and the relevant correlator in the 
chiral limit reads 

^ I d'xe^^^ {T3^^{x)3^{Q)) = {q,q, - q'g,,) n^+^(Q2) , = -q' . (53) 

Up to 0{al) the perturbative spectral function is given as 

F3 = 1.9857 - 0.1153n^, 

F4 = -6.6368 - 1.2001 0.0052 (54) 

where = 3 is the number of light quark flavors. The Wilson coefficient of the 
dimension-4 gluon condensate correction is known to order a^. 
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Since the correlator Il{Q'^)^~^^ is cutoff-dependent itself, we investigate the Adler 
function 

For the corresponding experimental V+A spectral function we have used the Aleph 
measurement |2S1 in the resonance region up to 2.2 GeV^. For the continuum region 
above 2.2 GeV^ we used 3-loop perturbation theory for as{Mz) = 0.118, and we have 
set the renormalization scale /i to Q. We note that the pion pole of the axial vector 
contribution has to be taken into account in order to yield a consistent description in 
terms of the OPE for asymptotically large Q. [2H] We have checked that the known 
perturbative contributions to the Adler function show good convergence properties. 
For Q between 1.0 and 2.0 GeV and setting = Q, the 3-loop (two-loop) corrections 
amount to 5% (7%) and 0.5% (1.4%), respectively. On the other hand, the 0{as) 
correction to the Wilson coefficient of the gluon condensate are between —16% and 
—6% for fi between 0.8 and 2.0 GeV. In analogy to the analysis for charmonium 
sum rules and setting Q = Q we have compared the local dimension 6 contributions 
contained in the running gluon condensate in Eq. ()55|) with the sum of all dimension 
6 terms in the OPE as they were determined in Ref. |28|. The latter consist of local 
four-quark condensates, which we have estimated by assuming vacuum saturation 
at the scale /i ~ 0.7 and using the parameters employed for the charmonium sum 
rules. We again found that the corresponding dimension 6 contributions have equal 
sign and roughly the same size. 

Our result for the running gluon condensate as a function of Q is shown in Fig.|3 
for 1 GeV< Q < 2 GeV. The shaded region represents the allowed range obtained 
from an analysis employing the full 0{al) perturbative contribution and the striped 
region is obtained at 0{a'j.). The uncertainties are due to the experimental errors 
in the spectral function and a variation of the renormalization scale /i in the range 
Q ± 0.25 GeV. The strong coupling has been fixed at as{Mz) = 0.118. We have 
restricted our analysis to the range 1 GeV < Q < 2 GeV because for Q ^ 1 GeV 
perturbation theory becomes unreliable and for Q > 2 GeV the experimentally 
unknown part of the spectral function at s > 2.2 GeV^ is being probed. The thick 
black line in Fig.0]shows the lattice-inspired running gluon condensate of Eq. (j52j) for 
Q = Q using the parameters from Sec. 15. II We find that our result for the running 
gluon condensate as a function of Q is consistent with a function that increases 
with Q. The result is also consistent with the lattice-inspired expression for the 
running gluon condensate, but the uncertainties of our extraction are too large and 
the running of the lattice-inspired condensate in the range between 1 and 2 GeV is 
too small to draw a more definite conclusion. 
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Figure 4: The running gluon condensate as a function of Q when extracted from 
the Adler function. The grey area represents the allowed region using perturbation 
theory at C(a^) and the striped region using perturbation theory at 0{a'j.). The 
thick solid line denotes the running gluon condensate from the lattice- inspired ansatz 
in Eqs. ior n = Q. 



6 Anomalous dimensions 

In the previous sections we have presented the formalism for the DOE assuming that 
the separation of fluctuations at different length scales can be carried out without 
ambiguities and that there is no need to introduce a regularization scheme. We have 
applied the DOE for QCD quantities, where nonperturbative effects are dominated 
by the gluon condensate, which can be defined such that its anomalous dimension 
vanishes to all orders. This allowed us to ignore anomalous dimensions at the level 
of precision intended in our investigations in Sec. El 

However, in general, the separation of fluctuations at different length scales is 
ambiguous and a regularization scheme is required to carry out the separation sys- 
tematically. In this general case, the delocalized condensates cannot be defined 
through the integrals in Eqs. (fT^ or (jTTj) because the resolution scale fl would itself 
serve as a UV-regulator and, therefore, interfere with the regularization scheme. 
Since the regularization of the short-distance coefficients in the DOE remains unaf- 
fected, this would lead to inconsistencies in the DOE. Thus, in the general case, the 
delocalized condensates are defined through the sum over local condensates given 
in Eq. (^01) (and its li-dimensional generalization), which should be considered more 
fundamental than Eqs. ()16|) and ^I7\i. We emphasize that with this definition there 
is, by construction, no interference between the DOE and the regularization scheme. 
Having in mind perturbative computations in quantum field theories in general, this 
means in particular that the delocalized expansion does - up to the summations that 
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define the short-distance coefficients and the matrix elements - not affect renormal- 
ization, i.e. the structure of anomalous dimensions, the running of coefficients and 
matrix elements, etc.. 

Coming back to the lattice-inspired expression for the running gluon condensate 
in Eq. ()52j) it is now clear that it was determined in the exponential cutoff scheme 
defined through Eq. ()27p . However, because the ^-dependence of the running gluon 
condensate is mainly due to the exponential behavior of g{x) and because the actual 
renormalization-scale-dependence of the running gluon condensate is small (due to 
the fact that it is generated by local condensates of dimension 6 and higher) the use 
of expression ()52j) for the comparisons in SecEl where the perturbative contributions 
were determined in the MS scheme, are justified. 

It is obvious that the delocalized expansion can in principle be applied to any 
quantum field theoretical computation, where a local expansion can be carried out. 
In particular, it can be applied in effective theories, and it allows for the systematic 
summation of matrix elements of operators with additional covariant derivatives in 
addition to the summation of large logarithmic perturbative corrections obtained 
from solving the renormalization group equations for the short-distance coefficients. 
In the following we briefly demonstrate the delocalized expansion in the presence of 
a regularization scheme - dimensional regularization to be specific - for fo/ fs at 
leading order in the heavy quark mass expansion. 

The meson decay constants are an important parameter in computations of mix- 
ing and lifetime-differences. They can be computed in heavy quark effective theory 
(HQET) in an expansion in AQCD/mq, where mq is the respective heavy quark mass. 
To be definite, we consider the decay constant of a pseudoscalar meson, which is 
defined as^ 

{Olqr^'^QlPiv)) = ifp^jM~pV^, (57) 

where q and Q are the light and heavy quark fields in full QCD, P{v) is the meson 
state with velocity v in full QCD, and Mp is the mesons mass. In the (local) l/mq- 
expansion the current q'j'^'-f^Q is written as a series of HQET currents of increasing 
dimension each of which is multiplied by a Wilson coefficient: 

q^^^Q ^ _ Co{mq, /i) q^^'^'^hq - Ci{mq, /i) q^^'^hq + .... (58) 

Here, hq denotes the heavy quark field in HQET and the ellipses represent contri- 
butions from currents with higher dimensions. The leading logarithmic expressions 
for the Wilson coefficients read [21] 

Co{mq,i^) = {^^y^^ C^{mq,^^) = 0, (59) 

where /Sq is the one-loop beta-function. Since the vacuum-to-meson matrix element 
of the current qj^j^hq is equal for B and D mesons at leading order in 1/mq, this 

^Since we discuss only the leading order in the l/mq expansion, our treatment also applies to 
vector mesons. 
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leads to the leading order relation 



fB 



Mb f asjmb) 



) 



_6_ 

25 



(60) 



for 4 active light flavors. Recent lattice measurements jSU] indicate that 



fB 



1.0 - 1.3, 



(61) 



but the RHS of Eq. gives /d/Zb ~ 1-5. It was found in Ref. [SI] that the dis- 
crepancy cannot be accounted for quantitatively by subleading l/niQ contributions, 
because the corresponding corrections are simply too large to allow for a quantita- 
tive improvement of the leading order relation in Eq. (j60p . In particular, for the D 
meson the local l/mg expansion seems to break down, and it appears that going 
beyond the leading order approximation for / fs is meaningless. 

We suspect that the same statement might also apply to the delocalized version 
of the l/mq expansion, but from the examinations in the previous parts of this 
paper one can expect that the delocalized expansion might provide a better leading 
order approximation for an appropriate choice of resolution scales. In the delocal- 
ized l/mq expansion the leading order short-distance coefficients are equivalent to 
the local ones shown in Eq. ()59|) . However, the vacuum-to- meson matrix element 
of the current Q'j^'j'^hQ becomes a resolution-scale-dependent quantity because it 
sums matrix elements of local currents of the generic form D^"g7^7^/iQ (n > 0). 
Evidently, the proper choice of the resolution scale is ~ mq. Thus the leading 
order expression for fn/fB in the delocalized expansion can be written as 



The common renormalization scale n in the f2-dependent matrix elements is set 
below the charm quark mass, but not fixed otherwise. Comparing expression (jU^ 
with the lattice result one finds p ~ 0.6-0.8. (Neglecting the anomalous dimension of 
the current q'j^'y'^hq gives essentially the same result due to the large uncertainty of 
the lattice measurements.) We are not aware of any other independent determination 
of the ratio p, but this scaling relation has an interesting physical interpretation 
in the framework of a model, where one assumes that the B and the D meson 
wave functions have approximately the same size^ and that the vacuum-to-meson 
matrix elements of the local currents D^^'q'y^'y'^hq scale like inverse powers of the 

^The results from the lattice simulations of Ref. [21] indicate that this assumption is not unre- 
alistic. This assumption is also consistent with results from potential model computations. |33| 




(62) 



where 



P 



(0|g7V K\D{v)){p,Q = m,) 
{OlqY'l" h\B{v)){p,n = nib) ' 



(63) 
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meson size. Assuming exemplarily that the meson wave functions have a Gaussian 
drop-off behavior (see e.g. Ref. |Mi) (0|g(-x/2)7V x/2) /iq(x/2)|P) - 
exp(— |xp/A^) one finds, using the 3-dimensional version of Eq. (|77j). 



1 + 



b 



1 + 



(64) 



A2r 



One can then derive that A ^ ~ 0.5-0.8 GeV, which gives an estimate for the meson 
size. On the other hand, for an exponential drop-off behavior ~ exp(— |x|/A) one 
obtains A'^ ~ 0.3-0.7 GeV. 



7 Summary and outlook 

In this paper we have proposed a new type of short- distance expansion of correlators 
of gauge invariant currents in QCD, which we name delocalized operator expansion 
(DOE). This expansion originates from nonlocal projections of gauge invariant cor- 
relation functions based on a delocalized version of the multipole expansion for the 
perturbatively calculable coefficient functions. Whereas the usual (local) OPE is 
based on a multipole expansion where the short-distance process in configuration 
space is written as the sum of delta-functions and derivatives of delta-functions, 
the DOE is based on a delocalized version of the multipole expansion using func- 
tions that have the width Q^^. We call fl the resolution parameter. In this paper 
we have constructed such a delocalized multipole expansion based on the Gaussian 
function and using the orthonormality and completeness properties of the Hermite 
polynomials. We emphasize that other, similar constructions are possible which 
have properties that are comparable to the one presented here and that they might 
be even better suited to certain problems than the one we have used here. 

In the DOE the condensates become quantities that are dependent on the resolu- 
tion parameter Q and are related to an infinite sum of VEV's of local operators with 
equal and larger number of covariant derivatives. The short- distance coefficients also 
become dependent on the resolution parameter Q and are related to a finite sum of 
short-distance coefficients to local operators with equal and fewer number of covari- 
ant derivatives. The relative weight of the local terms in these sums is governed by 
the resolution parameter Q. For Q ^ oo the DOE reduces to the usual local OPE. 
As is the local OPE, the DOE is parametrically an expansion in powers of A/Q, 
where A denotes the highest nonperturbative scale occurring in the corresponding 
gauge invariant correlation function and Q is the external short- distance scale at 
which the vacuum is being probed. If the resolution scale is chosen to be of order 
Q, there is an additional suppression of higher orders in the expansion by powers of 
a small number. This can be understood intuitively, because for f2 ~ Q the leading 
term in the delocalized multipole expansion represents a better approximation of the 
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short- distance process than the delta-function. In this way the DOE can account 
for non-local effects. 

Calculating the nonperturbative correction to the perturbatively determined 
ground state energy of a heavy quarkonium system, we have demonstrated the im- 
proved convergence properties of the delocalized expansion using a lattice-inspired 
toy model for the gluonic 2-point field strength correlator. Reversing the situation, 
we have extracted the running gluon condensate {{as/7T)G'^){fl) from experimental 
data using charmonium sum rules and the Adler function in the V + A channel of 
light quark pair production. We found strong evidence that the gluon condensate 
is indeed a resolution-dependent quantity and our results are consistent with recent 
lattice measurements of the correlation length of the gluon field strength correlator. 

The DOE can be applied in the framework of effective theories. By construction, 
the DOE does not interfere with the regularization scheme since the resolution- 
dependent short- distance coefficients and matrix elements are defined as sums of 
local short-distance coefficients and matrix elements, respectively. Thus, the DOE 
does not affect the renormalization properties of a theory. In the DOE we have 
derived the leading order expression for / fs in the heavy quark mass expansion. 

Further investigations and more detailed applications of the DOE are planned 
in forthcoming work. 



A Tensor structure of {g^GD'^G) 

Here we show that the lattice-implied dominance of the tensor structure in Eq. ()34|) 
for the gluon field strength correlator is consistent with a phenomenological deter- 
mination of the dimension 6 condensate {g'^G^^{0)DaD/sGf^x{0)) ■ When evaluated 
in Schwinger gauge x^A^ = |3] we are, according to Eq. concerned with the 
quantity l^g'^G'^^j^dpdpG'^f^ taken at x = 0. Following Ref. |21j it can be parametrized 
in Euclidean space-time as 

+ {^^lP^(xu + - 5f,aSu/3 " ^lyf^Spa) , (65) 

where 

= ^ (a'jX) ± ^ {9'fai.cGl^Gl,Gl^) , (66) 

and denotes a light flavor- singlet current. In Eq. (j65|) the tensor structure multi- 
plying 8 0~ is the same as the one resulting from an omission of the terms involving 
Di in the parametrization of Eq. (jHT|) . Now, we consider 

(g'G^dpdpG^) = 96 + 24 {Q- + 0+) , (67) 
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and we have to show that, phenomenologically, the contribution 9QQ is the domi- 
nant one, that is 

^°- 24(0-V) »^- 

Assuming exact vacuum saturation to express {g'^j'^) in terms of a square of the quark 

ya rib r^c \ 



condensate and using a, = 0.7 (/i = 0.7 GeV), ig^ f abcGl^Gl^Gl}i = 0.045 GeV^, 



and {qq) = —(0.24 GeV)^ as in Sec. 5.1, we obtain 



^4 AaAa 
y J aJ 



-2.3. (69) 



For the ratio Rq of Eq. this imphes 

Ro ~ 8.8, (70) 
which justifies the omission of the terms in the second and third hne of Eq. (j65j) . 

B Short-distance coefficients for quarkonium ground 
states energy levels 

f ( \ 702 m 

= 425 361^' 

^ , , 1108687904 

= 844421875 36¥ ' ^^'^ 
17608502062304825344 
^^^'^^ ~ 2270010244541015625 36 k^^ ' ^ ^ 

45707850430207738594209974779904 
~ 350274120613508115758056640625 36 A; ' ^ ^ 

3377025742042031485195560438171693564080529866752 
~ 754795366909616960779550333930811309814453125 36 ' ^ > 
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